Cold Nuclear Matter In Holographic QCD 



Moshe Rozali 1 , Hsien-Hang Shieh 1 , Mark Van Raamsdonk 1 , and Jackson Wu 2 

1 Department of Physics and Astronomy, University of British Columbia 
6224 Agricultural Road, Vancouver, B.C., V6T 1W9, Canada 

Theory group, TRIUMF, 4004 Wesbrook Mall, Vancouver, BC, V6T 2A3, Canada 



Abstract 

We study the Sakai-Sugimoto model of holographic QCD at zero temperature 
and finite chemical potential. We find that as the baryon chemical potential is 
increased above a critical value, there is a phase transition to a nuclear matter 
phase characterized by a condensate of instantons on the probe D-branes in the 
string theory dual. As a result of electrostatic interactions between the instan- 
tons, this condensate expands towards the UV when the chemical potential is 
increased, giving a holographic version of the expansion of the Fermi surface. We 
argue based on properties of instantons that the nuclear matter phase is neces- 
sarily inhomogeneous to arbitrarily high density. This suggests an explanation of 
the "chiral density wave" instability of the quark Fermi surface in large N c QCD 
at asymptotically large chemical potential. We study properties of the nuclear 
matter phase as a function of chemical potential beyond the transition and argue 
in particular that the model can be used to make a semi-quantitative prediction 
of the binding energy per nucleon for nuclear matter in ordinary QCD. 



1 Introduction and Summary 



QCD at finite temperature and chemical potential 

The phase diagram of QCD as a function of temperature and baryon chemical po- 
tential (or alternatively baryon density) displays a rich variety of phases and transitions 
(for reviews, see [HOE]). However, apart from the regimes of asymptotically large 
temperature or chemical potential, where some analytic calculations are possible, and 
of zero chemical potential, where reliable lattice simulations are possible, our knowl- 
edge of the phase diagram is based exclusively on extrapolations and semi-empirical 
toy models. For intermediate values of the chemical potential, numerical simulation is 
plagued by a notorious 'sign problem' (see for example [3]), while analytic calculations 
are not possible due to strong coupling. Thus, while there has been significant progress 
recently in understanding the qualitative features of the phase diagram, reliable quan- 
titative calculations that would definitively verify the proposed phase structure or de- 
termine the locations of various transitions or properties of the various phases seem a 
formidable challenge at present. A better understanding of the details of the phase dia- 
gram at intermediate chemical potential would have valuable applications, for example 
in understanding the physics of neutron-star interiors. 

Holographic models of QCD 

With the advent of the gauge theory / gravity duality [5], we have a new tool for 
studying the properties of certain strongly coupled gauge theories. While the original 
and most studied examples involve highly supersymmetric conformal gauge theories 
without fundamental matter, much progress has been made in constructing examples 
without supersymmetry [6], with confinement [6], with fundamental matter [7] and 
with chiral symmetry breaking [Hj. We now have examples of gauge theories with a 
known gravity dual that share most of the qualitative features of QCD, and the duality 
permits analytic calculations that would be otherwise impossible. 

It is obviously interesting to study these QCD-like theories in regimes for which 
neither analytic or numerical studies are currently possible in real QCD. One such 
regime is the near-equilibrium behavior of the theory at finite temperature . This has 
received a great deal of attention recently (see [I] for a review) since calculations in 
holographico models of QCD-like theories do a better job of explaining and predicting 
some properties of the quark-gluon plasmas produced in relativistic heavy-ion collisions 
than any other approach. In the present paper, our focus will be on another such regime 
as described above, the equilibrium properties at finite baryon chemical potential. 

There is already a large literature on studies of gauge theories at finite chemical 
potential using gravity duals (see [9] and references therein). Many of these consider a 
chemical potential for R-charge in theories with only adjoint matter. There have been 
some some studies of the behavior of theories with fundamental matter at finite baryon 
chemical potential, but the early examples of holographic theories with fundamental 

^^Here 'holographic' is a now conventional term referring to the equivalence between a higher- 
dimensional gravitational theory and a lower-dimensional field theory. 
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Figure 1: Type IIA string theory configuration for the Sakai-Sugimoto model. 



matter had both bosonic and fermionic fields carrying baryon charge. In these cases, 
the physics at finite chemical potential involves Bose condensation rather than the 
formation of a Fermi surface. In order to get behavior similar to real QCD, it is 
essential to study a theory with baryon charge carried exclusively by fermionic fields. 
Such a model was constructed a few years ago by Sakai and Sugimoto [ID] , and it is 
this model that we will focus on the present work. 

The Sakai-Sugimoto model 

The details of the Sakai-Sugimoto model are reviewed in section 2. Briefly, the 
model gives a holographic construction of a non-supersymmetric SU(N C ) gauge the- 
ory with Nf fundamental fermions. The gravity dual involves Nf D8-branes in the 
near-horizon geometry of N c D4-branes wrapped on a spatial circle with anti-period 
boundary conditions for the fermions. In the geometry, the compact direction of the 
field theory together with the radial direction form a cigar-type geometry, in which the 
D8-branes are embedded as shown in figure 1. The other directions include an S A car- 
rying N c units of D4-brane flux and the 3 + 1 directions of the field theory. In addition 
to Nf and N c , the theory has a dimensionless parameter A, the 't Hooft coupling at 
the field theory Kaluza-Klein scaled 

For small values of A, the scale Aqcd where the running coupling becomes large is 
well below the field theory Kaluza-Klein scale, and the low-energy physics should be 
precisely that of pure SU(N) Yang-Mills theory coupled to Nf massless (fermionic) 
quarkso Unfortunately, in this limit, the dual gravity background is highly curved so 
we are not in a position to study it. For large A on the other hand, the gravity back- 
ground is weakly curved, and so via classical calculations on the gravity side of the 
correspondence, it should be possible to map out the phase diagram of the field the- 

2 The model has another parameter, corresponding to the asymptotic separation between the D8- 
branes, but we focus exclusively on the case where the two stacks are on opposite sides of circle and 
extend down to the tip of the cigar. 

3 For recent work on adding quark masses, see [TTl 112) . 
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ory as a function of temperature and chemical potential and quantitatively determine 
properties of the various phases. 

We do not expect our results to agree quantitatively with real QCD (both because 
the Kaluza-Klein scale is not well separated from Aqcd for large A and because the 
classical calculations give only the leading terms in the 1/N expansion), but it would 
certainly be interesting to have a precise understanding of the phase diagram for a 
theory that is so similar to QCD. Indeed, at least some features of the phase structure 
and the qualitative behavior of certain transitions are likely to be the same as in 
QCD, and we might even hope for rough quantitative agreement for quantities that 
are relatively insensitive to A and N c (we will discuss one such quantity below) . 

The transition to nuclear matter 

Our focus in this paper will be on the part of the phase diagram for zero temper- 
ature and intermediate values of the baryon chemical potential. In real QCD, as we 
increase the chemical potential from zero, the equilibrium state (i.e. the ground state) 
continues to be the vacuum until some critical value of the chemical potential at which 
point it becomes advantageous for baryons to condense. A first approximation to this 
critical value is the baryon mass, since it is at this point where it becomes energetically 
favorable to add single baryons to the vacuum. In fact, the critical value is somewhat 
lower, since the baryons have a negative binding energy. At the critical value, we have a 
first order transition from the vacuum state to homogeneous nuclear matter with some 
minimal baryon density]^ The best estimate for the critical chemical potential comes 
by studying the masses of atomic nuclei as a function of nucleon numbers |15j . These 
are fit very well by the Weizsacker-Bethe semiempirical mass formula, which includes 
a term proportional to the number of nucleons, 

TTT'Vol = b vo lA 

to take into account the energy — b vo i due to strong interactions of each nucleon in the 
interior of a nucleus with its neighbors plus the average kinetic energy per nucleon 
(non-zero due to Fermi-Dirac statistics). The best fit for this energy is 

b vol = 16 MeV . (1) 

Ignoring electromagnetic interactions, this gives the binding energy per nucleon in 
the limit of large nuclei, and thus should be a good approximation to the value for 
homogeneous nuclear matter just beyond the transition. Thus, the critical chemical 
potential for the transition to nuclear matter in QCD should be approximately 

H C = M B (1 - ^-) « Mb(1 - 0.017) . 



4 It is important to note that we are talking only about QCD and ignoring electromagnetism here. 
With electromagnetic interactions, the binding energy per nucleon is actually greater in iron nuclei 
than in homogeneous nuclear matter, so the transition to nuclear matter is preceded by a transition 
to solid iron. 
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As we increase the chemical potential further, the baryon density and the energy per 
baryon will increase from their values just above the transition. Eventually we hit at 
least one more transition, to a phase characterized by quark-quark condensates pQ. 

In this paper, we will study the physics of the transition to nuclear matter in the 
Sakai-Sugimoto model at large A. Via classical calculations in the dual gravitational 
theory, we will be able to determine the critical chemical potential and calculate the 
baryon density and the energy per baryon ee(/z) for fi above the transition. 

Expectations at large N 

Since our gravity calculations will give results corresponding to the large N c limit 
of the field theory (with a fixed Nf), we should briefly recall the expectations for how 
baryons behave for large N c |13j . In this limit, baryon masses and baryon-baryon 
interaction energies go as N c , but the baryon size approaches a constant. Thus, we 
expect that both the baryon density above the transition and the binding energy per 
nucleon divided by the baryon mass to have a finite limit for large N c . These properties 
indeed follow from our calculations. 

One significant difference between the large N c theory and ordinary QCD is the ex- 
pected behavior at asymptotically large values of the chemical potential. In both cases, 
we have attractive interactions between excitations on the Fermi surface that result in 
an instability, but the nature of the resulting condensates is different. Whereas for 
N c = 3 the instability is a BCS-type instability, believed to lead to a color supercon- 
ductor phase, the dominant instability at large N c is toward the formation of "chiral 
density waves" [18], inhomogeneous perturbations in the chiral condensate with wave 
number of order twice the chemical potential. This suggests that the ground state for 
large N c QCD at large enough chemical potential is inhomogeneous, however the na- 
ture of the true ground state remains mysterious (see [TS] for a recent discussion). We 
believe that our analysis sheds some light on this question, as we will discuss shortly. 

Results for the Sakai-Sugimoto model 

In the Sakai-Sugimoto model, a chemical potential for baryon number corresponds 
to a nonzero asymptotic value of the electrostatic potential on the D8-branes, equal on 
both asymptotic regions of the D8-brane. Generally, this potential behaves asymptot- 
ically (for radial coordinate U to be described below) as 

A ~ fi B + E-^ + . . . . 

The baryon density ris is proportional to the asymptotic abelian electric flux E, so 
configurations with non-zero baryon density in the field theory correspond to D8-brane 
configurations with sources for the electric flux. These sources can be either string 
endpoints on the D8-branes which originate from D4-branes wrapped on the internal 
S A of the geometry [J4] or (for Nf > 1) configurations of the Yang-Mills field carrying 
instanton charge [161 EZJ- The latter can be thought of as the wrapped D4-branes 
dissolved into the D8-branes and expanding into smooth instanton configurations. 
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One flavor 

For any value of chemical potential, we always have a trivial solution for which 
the electrostatic potential is constant on the D8-branes and the baryon density is 
zero. However we can also consider translation invariant configurations with a uniform 
baryon density. In the single flavor case, which we consider first, the bulk description 
of baryons is in terms of pointlike instantons, since there are no large instanton config- 
urations in the abelian gauge theory of a single D8-brane. In this case, configurations 
with a uniform baryon density correspond to having some density of these pointlike 
instantons on the D8-brane. For a given value of the chemical potential greater than 
the critical value, we find some preferred distribution of charges on the D8-brane. The 
total baryon density for a given value of fj, may be read off from the asymptotic value 
of the electric flux, and the result increases smoothly from above the critical chemical 
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potential, approaching an asymptotic behavior n B oc //2. The charge distribution in 
the radial direction for a given value of \i represents the distribution of energies in the 
condensate of baryons in the field theory. In particular, the distribution has a sharp 
edge at some value of the radial coordinate which increases for increasing chemical 
potential, and this gives a bulk manifestation of the (quark) Fermi surface in the field 
theory. 

For the single flavor case, the transition to nuclear matter is continuous, unlike 
QCD, but it may be expected that the single flavor case is different due to the absence 
of pions which usually play a crucial role in interactions between nucleons. 

Two flavors 

In the case with Nf > 1, we can have nonsingular instantons on the Nf coincident 
D8-branes, and the minimum energy configurations for large enough /i are should in- 
volve smooth configurations of the nonabelian gauge field carrying an instanton density. 
While we might expect this to be homogeneous in the field theory directions, we argue 
that there are no allowed configurations of the D8-brane gauge field that are spatially 
homogeneous in the three field theory directions such that the net energy density and 
baryon density in the field theory are both finite. Thus, any phase with finite baryon 
density is necessarily spatially inhomogeneous. This has a simple interpretation: it 
suggests that at large N c , the nucleons retain their individual identities for any value 
of the chemical potential. Assuming that this holds true also for small A where the 
theory becomes 2 flavor QCD, this suggests that the chiral density wave instability of 
the quark Fermi surface in large N c QCD simply indicates that the quarks want to 
bind into nucleons even at asymptotically large densities. This is discussed further in 
section 5. 

To avoid the complication of directly studying inhomogeneous configuration, we 
approximate these by certain singular homogeneous configurations, arguing that our 
approximation should become exact in the limit of large densities. Within the context 
of this approximation, we study the behavior of the system as a function of chemical 
potential. 

Our model displays a first order transition to nuclear matter at some critical chem- 
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ical potential that depends on the parameter A, with the baryon density behaving as 
ns oc /i 3 for large \x. In the limit of large A, the critical value approaches the baryon 
mass, so the binding energy per nucleon is a vanishing fraction of the baryon mass at 
large A JU 

For large but finite A, we find the behavior 

/x c = M° (1 + ^ + ' 



A VAf 

where M B is the large A result for the baryon mass 

K = ^rM KK \N c . 
On the other hand, the baryon mass for large but finite A is [T6| [T7] 

M B = M2,(l + j + o(X 

It is interesting that the result for the binding energy per nucleon at the threshold for 
nuclear matter formation, 

E bind = M B - ^~M KK {d - c) , 

1 ITT 

is actually insensitive to the value of A for large A. Since we also know that this binding 
energy approaches some constant value in the limit of small A (the large N c QCD result 
with two massless flavors), then assuming a smooth behavior at intermediate values 
of A, we can treat the large A result as a prediction for the order of magnitude of the 
QCD result§ Noting that Mkk ~ J^qcd for large A, the value of the binding energy 
per nucleon extrapolated to N c = 3 becomes 

E bi nd = ^-Aqco(c' - c) « 7 MeV(c' - c) 

In order to reliably compute the the numerical coefficients c and c', we require knowl- 
edge of the nonabelian analogue of the Born-Infeld action, and (in the case of c') 
probably corrections to this involving derivatives of field strengths. However, assuming 
c' — cis of order one0we do obtain the same order of magnitude as the QCD result ([1]). 
We are not aware of any other methods to reliably estimate this binding energy from 
first principles, so it is possible that a more complete calculation in the Sakai-Sugimoto 
model would represent the most reliable analytic prediction of this quantity. 



5 While this statement is derived in the context of our approximation, we argue that it should be 
true in the full model. 

6 Another example with similar insensitivity to A for both large and small A is the free energy of 
N = 4 SUSY Yang-Mills theory. Here, it is indeed the case that the large A result for the free energy 
gives a good prediction of the order of magnitude of the the small A result (or vice versa) . 

7 We must also assume that our approximation scheme at least gets the right power of A in the 
correction to \i c . 
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Outline 



The remainder of the paper is organized as follows. In section 2, we review the 
Sakai-Sugimoto construction and collect various results necessary for our investigation. 
In section 3, we review the description of baryons in the Sakai-Sugimoto model and 
outline the basic approach for studying the theory at finite chemical potential. In 
section 4, we consider the single flavor case, calculating the baryon density as a function 
of chemical potential above the transition to nuclear matter. In section 5, we discuss 
the two flavor case, introduce our approximation, and set up a variational problem that 
determines the minimal energy configuration with a fixed baryon density (within our 
approximation). We then study the variational problem numerically for various values 
of chemical potential and baryon density to determine the critical chemical potential 
above which the minimum energy configuration has non-zero baryon density. 

Related Work 

Our work complements and extends various previous studies of the phase diagram 
for the Sakai-Sugimoto model. The behavior at finite temperature was analyzed in [22] • 
The behavior of the Sakai-Sugimoto model at finite chemical potential has also been 
discussed (with a different focus from the present paper) in [23l [2H [2E] . Discussions of 
the finite density behavior in other holographic models of QCD include [9| 1291 |28| [27 ] 134"] 

While this paper was in preparation, the paper [33] appeared, which has some 
overlap with the present work, in particular section 4.1. 



2 The Sakai-Sugimoto model 

The basic setup for the Sakai-Sugimoto model [10] begins with the low-energy de- 
coupling limit of N c D4-branes wrapped on a circle of length 2nR with anti-periodic 
boundary conditions for the fermions [6j. Apart from N c , this theory has a single 
dimensionless parameter 

the four-dimensional gauge coupling at the Kaluza-Klein scale. Because of the antiperi- 
odic boundary conditions, the adjoint fermions receive masses of order 1/R while the 
scalars get masses of order X/R due to one-loop effects. The coupling runs as we go to 
lower energies, becoming strong at a scale 

1 =* 

for some numerical constant c. As pointed out by Witten [5J, for small A, the dynamical 
scale Aqcd is far below the scale of the fermion and scalar masses and the Kaluza-Klein 
scale, so the dynamics should be exactly that of pure Yang-Mills theory. 

The field theory here is dual to type IIA string theory on the near-horizon geometry 
of the branes. The Lorentzian metric, dilaton, and four-form field strength are given 



7 



by 



3 . 3 



ds 2 = (jj^ 2 (f^dattfc" + f(U)dx 2 4 ) + (jf) 2 (j^dU 2 + U 2 dn 2 4 ) 



3 

u ^ 3 

e " = #s I 77" 



ft 



2tcN c 4 

F 4 — £4 

C^4 

where co>4 is the volume of a unit 4-sphere, e± is the volume form on S* 4 , and 

3 



f(U) = 1 - 



Uo 

u 



The X4 direction, corresponding to the Kaluza-Klein direction in the field theory, is 
taken to be periodic, with coordinate periodicity 2nR, however, it is important to note 
that this rr 4 circle is contractible in the bulk since the x 4 and U directions form a 
cigar-type geometry. 

The parameters i? 4 and U appearing in the supergravity solution are related to 
the string theory parameters by 



R\ = ng s N c l 3 s U = 



while the four-dimensional gauge coupling A is related to the string theory parameters 

as 

A = 2n? M 



R 

In terms of the field theory parameters, the dilaton and string-frame curvature at the 
tip of the cigar (the IR part of the geometry) are of order A 2 /N c and V% so as usual, 
supergravity will be a reliable tool for studying the infrared physics when both A and 
N c are large (in this case, with iV c >> A2). 

Note that this is opposite to the regime of A where we expect pure Yang-Mills theory 
at low energies. However, we may still learn about pure Yang-Mills theory by studying 
this regime, since many qualitative features of the theory remain the same and we 
might expect further that certain quantitative features may be relatively insensitive to 
the value of A (as for example with the free energy in M = 4 SYM theory). 



2.1 Adding fundamental matter 

Now that we have defined the adjoint sector of the theory, we would like to add fun- 
damental quarks. We keep the number of quark flavors fixed in the large iV c limit, but 
this means that the number of degrees of freedom in the fundamental fields (includ- 
ing the gauge field) is smaller than the number of degrees of freedom in the adjoint 
sector by a factor Nf/N c . Thus, for Nf fixed in the large N c limit, the influence of 
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the fundamental fields on the dynamics of the adjoint fields should be neglig ible@ In 
other words, what is known as the "quenched approximation" in QCD literature is 
exact in this limit. This implies that adding the additional matter does not modify 
the geometry, and indeed the construction of Sakai and Sugimoto (following earlier 
constructions) involves adding branes to the geometry which are treated in the probe 
approximation. 

The Sakai-Sugimoto construction is motivated by the observation that the light 
open string modes living at a 3+1 dimensional intersection of D4-branes and D8-branes 
give rise to chiral fermion fields on the intersection without accompanying bosons. 
Thus, to the original D4-branes, which we can take to lie in the 01234 directions with 
the £4 direction periodic, Sakai and Sugimoto consider adding a stack of Nf D8-branes 
and a stack of Nf anti-D8 branes separated at fixed locations in the X4 directions and 
extended along the remaining directions. This configuration is unstable before taking 
a near horizon limifl nevertheless, one can obtain a stable configuration of the probe 
branes in the bulk geometry by fixing the asymptotic positions of the D8 and D8-bar 
stacks in the X4 direction. The X4 positions of the branes are free to vary as a function 
of the radial direction U in the bulk of the geometry, and charge conservation implies 
that the two stacks necessarily join up in the interior of the geometry. Thus, (in the 
zero-temperature situation that we are considering) we really have just a single set of 
D8-branes, bent so that the orientation in the two asymptotic regions is opposite (see 
figure 1). 

The specific embedding of the D8-branes in the bulk depends on the asymptotic 
separation of the stacks (and also any distribution of matter on the branes), but we 
will focus exclusively on the case where the two asymptotic parts of the D8-brane stack 
sit at opposite sides of the D8 circle, in which case each side simply extends to the tip 
of the cigar along a line of constant X4 as shown in figure 1. The corresponding field 
theory has all flavors massless. 

2.2 D8-brane action 

To understand the physics of the probe D8-branes, we will need the action for the 
worldvolume D8-brane fields in the background above. We will begin by discussing the 
action for a single D8-brane before discussing the nonabelian generalization. 

The Born-Infeld action for the worldvolume D8-brane fields (in the case of a single 
brane) is 

S = -/i 8 J &<TfT*\f- det(g ab + F ab ) 

8 It would be quite reasonable to argue that we should keep Nf/N c fixed for large N c to obtain a 
theory that is most qualitatively similar to QCD, since then the number of degrees of freedom in the 
adjoint and fundamental sectors of the theory remain of the same order of magnitude for large N. 
However, this limit is much more difficult to study using supergravity, since then the back-reaction of 
the matter branes, to be described presently, must be taken into account. 

9 This instability is actually absent in the case we consider the stacks sit at opposite sides of the 
circle 
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where 

F = 2na'F 

We also have a Wess-Zumino term 

Here, only the C 3 term contributes. Noting that F 3 is the derivative of the five- 
dimensional Chern-Simons form, u$ and integrating by parts, we get 

S = -/i 8 / F 4 A co 5 . 



After integrating over the sphere, this gives 



N 

S = -JL / UB (A) (2) 



24tt 2 

where du^ = F A F A F. For a single D8-brane, W5 = i A F A F. 

To simplify the Born-Infeld action, we can choose to identify the worldvolume and 
spacetime coordinates in the sphere and the field theory directions, and parameterize 
the profile of the brane in the U and X4 directions by U(cr) and X(a) respectively (we 
will soon focus on the solution where X(a) is constant). 

We will be interested only in time-independent configurations homogeneous and 
isotropic in the spatial directions of the field theory (which we label by indices k). 
The most general configurations we will consider will have non-zero F ai , Fij, and F 0a , 
all functions only of a. 

Integrating the determinant from the sphere directions over the sphere, we get a 
factor 
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while the remaining five- dimensional determinant is 

- det(g^ + F^) = -(G 00 g aa + F 2 CT + g 00 F at (g + F) ij F aj ) det(G^- + F i3 ) 

with 

g aa = Gud a Xd a X + Guud^Ud^U . 

Note that we are using Gjj here to refer to the spacetime metric and g ab for the 
worldvolume metric. The final result (in the Abelian case) is 

S DBI = -^- n 2 Rl I d\vdaV 
9s 3 





F ai + (l^ijkFiaFjk) 2 \ (3) 

This action is manifestly invariant under reparametrizations of a. The nonabelian gen- 
eralization of this action is known only up to F 6 terms. Up to order F A , we symmetrize 
all of the nonabelian field strengths in expanding the square root and take an overall 
trace. However, this symmetrized trace prescription is known to fail beyond order F A . 
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2.3 Chemical potential for baryon charge 

We would like to study the theory at finite chemical potential for baryon charge or 
alternatively, the theory with a modified Hamiltonian density 



H = H + fxB 



where B is the baryon charge density operator 

B = B L + B R = ^y L + ^ R . 

This is equivalent to adding a term —fiB to the action since there are no time derivatives 
in B. Turning on the operator B in the boundary gauge theory with real coefficient 
H should correspond to turning on some (real) non-normalizible mode in the gravity 
picture. From the original brane setup, we know that the operators Bl and Br couple 
to the time-components of the D8 and D8 brane gauge fields respectively. We will 
see below that the equations of motion for these fields require them to approach some 
constant values in the UV part of the geometry. If we describe the probe branes as 
above with a single gauge field for the whole configuration, then we have two such 
constant values, 

Aoo = A (a = oo) 

and 

A_oo = A (a = -oo) 

These two values give the chemical potentials for the operators Bl and -Br@ Thus, to 
work at finite chemical potential for baryon number, we require that the value of A 
in both asymptotic regions of the D8-brane approaches the constant lib- 



2.4 Asymptotic solutions 



In the simple case where the D8-brane is at constant X4 and we assume that only the 
electrostatic potential is turned on, the Born-Infeld action above reduces to 



S 



DBI 



9s 3 4 



dad A xU^ 



1 



lf(U) ° 



BJJdJJ - d„Ad n A 



(4) 



The reparametrization invariance allows us to chose U(a) to be whatever we like. For 
a given choice of U, the equation of motion for A away from any sources (which we 
assume are localized in the infrared part of the geometry) is 



0. 1 ^r\u\ 

9s 3 



f(U) a 



dJJdJJ - d^Ad^A 



d n A 



(5) 



10 We give an argument in appendix A to establish that Bl and Br are turned on with the same 
sign if Aoo and A-oo have the same sign. 
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The quantity in round brackets is analogous to the conserved electric flux. Integrating 
and rearranging, and choosing o = U (valid for either half of the brane), we get 

E 

d u A = = , (6) 

y/f(U)(U* + E?) 

where E is an integration constant proportional to the conserved flux. Solving this, we 
find 

r-oo 

A = Anr, — I du- 



ll 



a/F(m)(m 5 + E 2 ) 



A 2 E 

— Aao + -— J + • • • 

6 {J 2 

valid in the region outside the sources. The constant E is the normalizible mode of A 
in the asymptotic solution, so the values of E for the two sides of the brane correspond 
to the expectation values for Bl and Br in the field theory. 

3 

In general, the sum of the Es for the two halves of the brane (times y-^n 2 Rl (2na')) 
is equal to the total charge density on the brane, 



^8 8 
9b 3 



7i 2 Rl{2na'){E 2 + E 1 ) = q 



If we fix A^ = A_ OQ as we have argued corresponds to a chemical potential for baryon 
number, and we assume that the sources are symmetric under a reflection in the a 
direction, then for continuous A we must have Ei = E 2 , and 

^7r 2 4(2W)£ = g/2 (7) 
9s 3 

Since the charge density in the bulk (divided by iV c ) corresponds to the baryon density 
in the field theory, we obtain 



3 Baryons 

We have seen that configurations with non-zero baryon charge density (as measured 
by the asymptotic electric flux E) require sources for Aq on the D8-branes. The basic 
source for Aq is the endpoint of a fundamental string. In order to have some net charge, 
we need the number of string endpoints of one orientation to be unequal to the number 
of string endpoints of the other orientation. So we need a source for fundamental strings 
in the bulk. In our background, such a source is provided by D4-branes wrapped on S 4 
|14j . These necessarily have iV c string endpoints, since the background D4-brane flux 
gives rise to N c units of charge on the spherical D4-branes, so we need N c units of the 
opposite charge (coming from the string endpoints) to satisfy the Gauss law constraint. 
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Thus, we can get a density of charge on the D8-brane by having a density of D4-branes 
wrapped on S 4 in the bulk, with N c strings stretching between each D4-brane and the 
D8-brane. 

In the case where we have Nf > 1 D8-branes, there is another possible picture of 
the configurations with baryons [TBI E] • To see this, note that a D4-brane / D8-brane 
system with four common worldvolume directions is T-dual to a D0-D4 system. In that 
case, it is well known that the DO-branes can "dissolve" in the D4-branes, where they 
show up as instanton configurations of the spatial non-abelian gauge field. Similarly, 
our baryon branes can dissolve in the D8-branes (if we have Nf > 1) and show up as 
instantons. Indeed, the Chern-Simons term (T5]) gives rise to a coupling 

S = ^ J A Tr(F A F) (9) 

between the instanton charge density and the abelian part of the gauge field, showing 
that instantons act as a source for the electrostatic potential on the branes. 

The question of which of these two pictures is more appropriate is a dynamical one, 
but it turns out that the dissolved instantons give rise to a lower energy configuration 
since the electrostatic forces prefer the instanton density to be delocalized [TBI [T7]. 



3.1 Baryon mass 

The baryon mass was estimated originally by Sakai and Sugimoto [1 Oj as the energy of 
a D4-brane wrapped on S* 4 and located at the tip of the cigar. Since we will also need 
to know the potential energy for such branes, we briefly recall the calculation. Starting 
with the Born-Infeld action for a D4-brane wrapping S A , 



S = -fi 4 d^e-W-det(g ab ) 



and integrating over the sphere, we get 

Sda = --l* 2 Rl [ dtU(t) (10) 

g s 3 J 

as the velocity independent term in the action (the negative of the potential energy). 
The minimum energy occurs for U = Uq, and this gives the baryon mass 

M ° B = ^-7T 2 RlU = —-XN C 
B g s 3 4 27nR 

This agrees with the Yang-Mills action for a pointlike instanton configuration on the 
D8-brane [16] . Both of these calculations ignore the energy from the electric flux 
sourced either by the string endpoints coming from the wrapped D4-brane or by the 
instanton density. To take this into account, the authors of [16j and p2] considered 
more general smooth instanton configurations with varying scale factor, inserting these 
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into the Yang-Mills approximation to the D8-brane action. They found that the optimal 
size for the instanton behaves as A~2 ; and that the baryon mass is 



This method ignores the effects of the non-trivial geometry on the Yang-Mills con- 
figuration and also does not include effects from the a' corrections to the D8-brane 
effective action, which should be important, since for large A, the instanton is small so 
that derivatives of the Yang-Mills field strength are large. Thus, as the authors point 
out, the numerical coefficient d should probably not be trusted. On the other hand, an 
analysis of the effects of Born-Infeld corrections [Hj indicates that at least the power 
of A in the correction to the mass and in the instanton size should be reliable. 

3.2 Critical Chemical Potential 

We have seen that turning on a chemical potential in the gauge theory corresponds to 
including boundary conditions Aq = fi for the two asymptotic regions of the D8-brane. 
For any /i, one solution consistent with these boundary conditions is to have constant 
A Q everywhere on the brane. This represents the vacuum configuration in the field 
theory. However, beyond a certain critical chemical potential, this solution is unstable 
to the condensation of baryons. 

The critical value of the chemical potential should not be larger than the baryon 
mass. At this value, a zero-momentum baryon has effectively negative energy in the 
modified hamiltonian, so it is advantageous to add baryons to the vacuum. If there 
were no interactions between the baryons, the critical chemical potential would be 
exactly the baryon mass. Note that even in the absence of interactions, the baryon 
density above the transition is limited by the Fermi statistics for the baryons for odd N 
or in any case by the Fermi statistics of the quarks. The condensate will have occupied 
all states whose Fermi energy is less than the chemical potential. In this case, the 
baryon density will rise smoothly from zero above the critical chemical potential and 
the transition will be second order. 

With short range repulsive interactions, the story would be qualitatively similar, 
with a slower growth in the baryon density as the chemical potential is increased. In 
QCD, however, we have attractive interactions, and this lowers the critical chemical 
potential below the baryon mass. With the repulsive interactions, there is a specific 
nonzero value of the baryon density for which the energy per baryon is lowest, and 
when the chemical potential is increased to this value the baryon density jumps from 
zero to this density. 

In the next sections, we will study this transition to nuclear matter in the Sakai- 
Sugimoto model for one flavor (section 4) and two flavors (section 5). In the first case, 
it appears that the transition is second order, unlike QCD, while in the multi-flavor 
case, we find some evidence for a more realistic first-order transition. 
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4 One flavor physics 



In this section, we study the physics of the Sakai-Sugimoto model at finite chemical 
potential in the simpler case of a single quark flavor. Here, we have only a single 
D8-brane in the bulk, and we can use the abelian Born-Infeld action for our analysis. 
Since the abelian gauge theory does not support large instantons, the wrapped D4- 
branes cannot dissolve into the D8-branes, so the baryons are pointlike charges on the 
D8-brane that source the electrostatic potential. For chemical potential larger than the 
baryon mass, it is favorable for some of these baryons to condense, and we would now 
like to determine the baryon density as a function of chemical potential for /i above 
the critical value. 



4.1 Localized source approximation 

As a first approximation, we make the simplifying assumption that all the pointlike 
instantons sit at U = Uq. More realistically, the charge should spread out dynamically, 
via electrostatic repulsion; we will include this effect in section 4.2. 

In our simple approximation, the relevant action is the Abelian Born-Infeld ac- 
tion (j3J), together with the action taking into account the baryon masses and their 
interaction with the electromagnetic field on the brane. 

i 



jjg 8 



3 



S = -—-n z Ri / dUdrxU* 



9s 3 



1 d~Ad a A 



f(U) 
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where the terms in the last line are the potential terms taking into account the en- 
ergy from the charges in the electrostatic potential and the masses of the pointlike 
instantons. 

To obtain the energy, we perform a Legendre transform, but it is convenient first 
to rewrite the first term in the second line as 

^A(U ) = - ( dU^A'(U) 

2na' V ; 2na' J 2na' V ; 

since we will be holding A(oo) = fi fixed. Performing the Legendre transform (which 
amounts to taking the negative of the action, since we are only looking at static con- 
figurations), and rewriting everything in terms of the electric flux (J6]), we find 



Zfiux = 2 ■ — ^ 2 Rl [ dll^-=(\ l + ^-r - 1) - (//-// c )n B 
9s % Juo Vf V U 5 



If 7 

——TT 2 RjU^h(e) -(ji- fi c )n B 
9s <j 



5 

where we have defined e = E/Uq and 

/oo 
dx(y x 5 + e 2 — x l 



V 1 - V* 3 
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In the first term, we have included a factor of 2 to take into account the energy f rom 
both halves of the D8-brane. 

For p > p c , the combined energy from the string endpoints (or Chern-Simons 
action) and the D4-brane mass (or Born-Infeld energy of the instantons) is negative 
and should be proportional to ns, while the energy from the flux is a positive function 

7 

of ub which behaves as n 2 B for small n# and n B for large n#. Thus, there will be some 
positive value of ns where the total energy is minimized. 
Defining 

6na'p 

so that p = 1 corresponds to p = M B , and using the relation ([7]) between n B and E, 
the total energy may be written as 



S = J s Y n2R * U ° (Me)-^(A-l)e) 



From this, we find that the energy is minimized when 

1 -(p-l) = h\e). 

This can be inverted to determine the relationship between ns (proportional to e) and 
p above the transition. For small p — p c , we find 

e ~ — (p, — 1) 

7T 

SO 

ns oc /i — fi c small \i — \i c . 

For large \x we have 

e ~ 0.021/23 

so 

n B oc /*2 large \i — \i c 



4.2 Dynamical charge distribution 

The analysis of the previous section assumed that all charges were localized at U — Uq. 
Presumably, the charges would prefer to spread out dynamically. To take this into 
account, we can define a charge distribution Pb(U) which we would like to determine. 
For a given p, the action is given in terms of a Lagrangian density 

c = ~ cul (ih) ~ d ' Ad - A ) 1 + ^ ApB ~^ UPB - 

where 

c = ^ S -M . 

9s 3 
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Here, the second term is the action arising from the string endpoints, while the third 
term takes into account the potential energy from the baryon masses (recalling that 
the action for a wrapped D4-brane at location U is proportional to U). 

For a given ps, the electric flux is determined by solving the equation of motion for 

A, 

Vs8 - a 



{2ira')du | ^ 2 R\m 
9s 3 



' - duAduA 



f(U) 



2 



duA = p B (U)N c (11) 



This gives 

Pb(U) = K J duE 
where we have defined an electric flux 

E(U) = ul(j±--(d u Af)^d u A. 

We can now reexpress all terms in the action in terms of E and Legendre transform 
(which again amounts to switching the sign) to find the energy. We obtain 



£ 

2C 



dU 

U 



-^{x/U^ + E 2 -1/3) + -UduE 

V j 3 



where we have included an extra factor of 2 in the denominator on the left side since 
we are integrating over only half the brane on the right side. To maximize this, we can 
first minimize over all E(U) such that E{Uq) = 0, E(U — > oo) = E^, and djjE > to 
determine £(E 00 ,fi). Then we can minimize over E^. 

Varying the energy functional with respect to E, we find that the energy functional 
is locally stationary if and only if 



e _ (12) 



(U 5 + E 2 )i 3 

This satisfies E = for U = (7o as desired but approaches arbitrarily large values for 
large U. On the other hand, our constraints dyE > and E(U — > oo) — > Eoo imply 
that E can never exceed E^. It is straightforward to check that the local contribution 
to the energy from a point U is a function of E that decreases from E = to the 
optimal value (fT2l) and then increases again, so when the value f|T2|) exceeds Eoo, the 
best we can do to minimize the energy is to set E = E^. We conclude that the 
minimum energy configuration for fixed p and fixed E^ is 

5 

E = /- U < U max 

ypi max (13) 

E = Eoo U > U m ax 

Here U max represents the extent of the charge distribution in the radial direction, and 
is related to E^ as 

-^oo _ \/ f(U ma x) 



(U 5 4 E 2 )3 

max ' oo/ 



(14) 
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We can now write the energy as a function of Eoo, or more conveniently, Uoo as follows. 
We define a function g(x) by 

5 

X 2 

9 j 



/(*) 

where 

/» = 1 - \ > 

and define 

= (\/a: 5 + t/ 2 . 

v/( x ) 

Then in terms of u — U max /Uo and , the energy is given by 

i ( f u f°° 1 f u 1 1 

£ = 2CU$ | J dxH(x, g{x)) + J H(x, g(u)) - -J g(x)dx + -ug{u) - -p,g(u) 

where as in the previous section, we define 

We can now minimize this as a function of u. The result is 



poo 

ji = u + 3 / dxd g H(x,g(u)) 

J u 



To compare with the results of the previous section, we note that (using (fill) ) the 
dimensionless variable e proportional to the baryon mass is related to u by 



/(«) 



[u 5 + e 2 )2 -3 
From these, we find that for small u — 1, 

/} — 1 = ci(it — 1)5 c x « 1.814 small it — 1 



or 

1 . 

e ~ — (u — 1_) 

7T 



where we have used (114ft . Thus, for small e we obtain the same result as in the previous 
approximation, with 

n B oc(iJL- fi c ) 

for small \i — /i c , where the critical value of \x is as before. For large /i, we find 

\i — > c 2 u c 2 ~ 1.697 large u 
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or 

e ~ 0.0942/tf . 

Again, we find that 

tib oc /i 5 . 

Thus, the qualitative behavior of is the same as in the simplified model of the 

previous section, though the numerical coefficients come out different. We also found 
the behavior of the energy density: 

£ oc (/i - yu c ) 4 

for (/i — /i c ) small, and 

8 oc M 7 / 2 

when n is large. 

It is interesting that (in this approximation) the charge distribution has a sharp 
edge at U = U max which progresses further and further towards the UV in the radial 
directions as the chemical potential is increased. In the field theory picture, the radial 
direction represents an energy scale, so the charge distribution we find in the bulk 
should be related to the spectrum of energies for the condensed baryons. The edge of 
the distribution is then a bulk manifestation of the Fermi surface. 

Since our large N c calculation does not distinguish between even and odd values of 
N c , it is insensitive to whether or not the baryons are fermions or bosons. Thus, the 
Fermi surface that we see should probably be thought of as the quark Fermi surface. 
It is interesting that the fermionic nature of the quarks in the field theory arises in the 
bulk from the classical electrostatic repulsion between the instantons. 

5 Two massless flavors 

For Nf = 2, the authors of [T6J [17] argued that single instantons on the D-brane prefer 
to grow to some finite size on the baryon in order to balance the electrostatic forces 
which tend to make the instanton spread out with the gravitational forces which prefer 
the instanton to be localized as much as possible near the IR tip of the D8-branes. 
From these considerations, we also expect that the minimum energy configurations 
with nonzero baryon density will involve some smooth configuration of the nonabelian 
gauge field on the D8-brane locally carrying an instanton density Tr(F A F). In this 
section, we consider such configurations. 
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The absence of homogeneous configurations 

We first consider static, spatially homogeneous configurations, such that is trans- 
lation invariant in the 3+1 directions of the field theory and rotationally invariant (up 
to a gauge transformation) in the three spatial directions (which we denote by an index 
i). The general configuration of the spatial gauge field with these symmetries is 

A a = A i = -L- ai h{a) (15) 

47r<x 

for an arbitrary function h(a). These givd3 

Fij = --^—e ijk a k h 2 (o) F ia = --aiti(a) . (16) 

47r<x 2 

From these, we find that 

We see that unless both h and h! vanish for o — > ±oo, the Yang-Mills action density 
integrated over a will diverge, corresponding to an infinite energy density in the field 
theory. On the other hand, we find 

(F A F) 123<7 = -^-h\o)h>{a) = ^(liV))l JxJ • 

In order that we have a configuration with finite baryon density in the field theory, 
we require that this instanton density, integrated over the sigma direction be non- 
zerc0. But this requires that h(oo) ^ h(—oo), and we have already seen that such a 
configuration will result in an infinite energy density in the field theory. 

The apparent conclusion for the dual field theory is that there are no spatially ho- 
mogeneous configurations with finite non-zero baryon density and finite energy density. 
Now, there certainly are non-homogeneous configurations with finite average energy 
density and finite average baryon density: we can simply take a periodic array of indi- 
vidual instantons. For large enough chemical potential (greater than the energy density 
divided by the baryon density), such configurations are favored over the vacuum, so 
we will certainly have a phase transition to a phase with nonzero baryon density as 
the chemical potential is increased. However, our observation suggest that this phase 
cannot be spatially homogeneous. 

Interpretation of the inhomogeneity and origin of the chiral density wave 

The inhomogeneity of nuclear matter is not unexpected, and indeed is what we 
have for real nuclear matter at low densities (e.g. in the interior of large nuclei). 

11 We use conventions where {<jj, <Tj} = 25^-1 and recall that F = (2na')F. 

12 To see this, note that the abelian electrostatic potential Aq couples to Tr(F A F), so that the 
change in the action upon a constant shift in Aq (corresponding to a change in the baryon chemical 
potential) is f daTr(F A F). 
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It simply reflects the fact that the individual nucleons retain their identities (and 
therefore that the baryon density is clumpecQ). What is perhaps surprising is that 
the inhomogeneity seems to have a topological rather than a dynamical origin from the 
bulk point of view, following from basic properties of instantons. It follows that even at 
arbitrarily high densities, the nuclear matter will be inhomogeneous, though the scale 
of the inhomogeneities should become shorter and shorter as the instantons pack closer 
and closer together. This suggests an interpretation of the DGR "chiral density wave" 
instability of the quark Fermi surface [18] at asymptotically large chemical potential: 
that even at arbitrarily high densities, quarks in large N c QCD bind into distinct 
nucleons, in contrast to the quark matter phase with homogeneous condensates that 
we expect at large /i for finite N c . This may be related to the property that the density 
of a baryon diverges for large N c and thus the baryon is more and more sharply defined 
in this limit. 

Our approximation 

The absence of homogeneous configurations with finite baryon density complicates 
the analysis of the phase transition and the properties of the nuclear matter phase. We 
will not attempt to study the inhomogeneous configurations directly here. Rather, we 
will describe an approach that approximates the inhomogeneous configurations with 
singular homogeneous configurations. 

Our approach is motivated by the observation that in the limit of infinite baryon 
density, the bulk configuration should become homogeneous. Such homogeneous con- 
figurations are singular at the core, corresponding to a divergence of the instanton 
charge density. For example, we can have a self-dual configuration of the form (fT5|) if 
we choose 

%) = - • (17) 
a 

This should arise from the limit of a periodic array of instantons for which the separa- 
tion is taken to zero while adjusting the scale factors to yield a non-trivial configuration 
in the limit. We expect that some similar configuration^! should arise in our case as 
the minimum energy configuration in the limit of infinite chemical potential. 

As we move away from infinite density, the minimum energy configuration will only 
be approximately homogeneous. We expect, however, that the averaged field strengths 
and instanton density should be qualitatively similar to those for the configuration (fTTl) 
but with finite values at a = 0. This behavior can be achieved in a configuration of the 
form ffTBT) for which h is an odd function like f[T7|) but with some finite limit at a = 0. 
Such configurations are singular at cr = 0, but we will ignore any effects associated 
with the singularity at a = since we are using our configurations to approximate non- 
singular inhomogeneous configurations that do not have any pathologies at a = 



13 Quantum mechanically this would be reflected in the behavior of density-density correlation func- 
tions. 

14 not necessarily self-dual since we are working with the D-brane effective action in a nontrivial 
geometry 

15 This is similar in spirit to replacing a nonsingular charge distribution with a localized singular 
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In particular, we might expect that our approximation becomes exact in the limit of 
infinite baryon density where we can have homogeneous configurations. We will find 
evidence below that supports the validity of this claim. More generally, we find results 
that are in accord with various physical expectations, providing further evidence for 
usefulness of our approximation. 



5.1 Energy density for approximate configurations 

We would now like to analyze the behavior of the model as a function of chemical 
potential in the approximation where we consider only configurations of the form (|15|) . 
taking h to be a monotonically increasing function for a > that takes some finite 
(negative) value at a = and vanishes for a — > oo. In practice, we work with the action 
for half the brane, assuming that h is an odd function so that all the field strengths 
are symmetric about a = 0. As we mentioned above, such configurations are singular 
at y — but we ignore any effects of the singularity, motivated by the expectation 
that the nonsingular contributions may provide a good approximation to the averaged 
quantities for the non-singular inhomogeneous configuration that we should really be 
studying. 

The configuration of the spatial SU(2) Yang-Mills field carries instanton density, 
and therefore acts as a source for the abelian electrostatic potential on the D8-branes. 
In order to determine the potential A(U) for a given h(U), we need the equation of 
motion for A, which should come from the non-abelian generalization of the Born-Infeld 
action ([3]) and the Cher n- Simons action (Q. 

As we have noted, the nonabelian generalization of the Born-Infeld action ([3]) is 
known only up to F 6 terms. In the absence of the full result, we will work with a naive 
ordering prescription in which we simply insert our ansatz into the abelian expression 
(J3J) and (noting that each product of Fs above gives an identity matrix) evaluate the 
trace. This will give us results that are precisely correct in the limit where the field 
strengths are small and only the Yang-Mills terms in the action are important, but we 
should not trust numerical coefficients whose calculation depends on the higher order 
terms in the Born-Infeld action. 

Inserting the ansatz ( Tl6i) into ([3]), we find (in the o = U coordinates): 



4(2tt«') 2 
(18) 



while the Chern-Simons term (j2j) gives: 



N 

S = — c - I Tr(A A F A F) 



24vr 2 
N, 



c 



1287r 6 (a')' 



dUAdu(h 3 (U)) . (19) 



distribution with the same multipole moments. 
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If we define 
and 



F = U x l{U/R A f J- J^Sl 



(Ana'] 



2 



then the action takes the form 



S = -C J dUFsjG - (duA) 2 + kj Adu(h 3 ) 



where 



k 



N, 



c 



1287r 6 (a') 4 

and 



C = l X 2 >^R 



3 g 8 



The equations of motion for the electrostatic potential A are 

Cd v E = kd v {h 3 ) 



where 



From this, we conclude that 



FduA , . 

E = = . (20) 

^G-(duAf 



kh 3 = C{E - E^) (21) 



where we have determined the integration constant by demanding that h vanish as 
U — > oo, as is required for finite energy configurations. Since E vanishes by symmetry 
at U = Uq (assuming that there is no delta function charge distribution at U — U ) we 
see that the asymptotic value of E is related directly to the value of h at U = Uq by 



khl = -CE^ . (22) 

We may therefore rewrite (I2ip as 

E = ^{h 3 -hl) 

Using this result, the electrostatic potential may be determined in terms of h by in- 
verting (1201) . 

We may now write an expression for the energy density of a configuration for a 
given value of h(U). 
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Starting with the actions ([TBI) and (|T9l) . we can derive the 3+1 dimensional energy 
density via a Legendre transformation as we did earlier. We find 



e = c du 



FJG-{duAf-F h=0 ^G^ 



+ k / duA(h 3 - hi) - kA^hl 



where we have subtracted off the energy density of the unexcited brane such that the 
vacuum state is normalized to zero energy. We can now rewrite the energy in terms of 
h, assuming that the equation of motion for A is obeyed. We have first 

S = C U dU(^G(F* + EZ) - F y/G^\ - CA^E^ 



Now writing E in terms of h as above, changing variables to x — U /Uq, defining 

h 



(23) 



An 



and 



we finally have 



An \i 

3 M^=°° 



(24) 



£ 



CU3 



Rt 



dx 



/(*) 



5 

12 



(25) 



Using the definition (1231) . and the relations (1221) and ([7]), we find that yo is related to 
the baryon density by 



n B 



7T / 4 g s NJ s 



9tt R V ° 



12^3 V 

- -?-aV- 

~ 27^ AoVo R s 

Thus, minimizing this expression for jj, = and fixed yo will give the minimum energy 
density for a fixed baryon density, which we denote by 



£min(yo) 

The energy density per baryon is then proportional to £ m inlv\i an d as we have argued 
above, the minimum of this gives the critical chemical potential. In the next section, 
we will analyze the functional (125]) . to obtain results for the behavior of £ m i n (yo) and 
for the critical chemical potential. 
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5.2 Results 



In this section, we discuss the evaluation of the baryon density for a given chemical 
potential based on minimizing the energy functional (1251) . Demanding that the func- 
tional is stationary under local variations of y gives a second order differential equation 
for y. For a given initial value y^ we find that there is a particular value of the ini- 
tial slope y' for which the solution approaches as x — > oo. For larger or smaller y' Q 
the solution approaches positive or negative infinity respectively for x —>■ oo, giving a 
diverging energy functional, so the minimum energy configuration must correspond to 
the solution with boundary condition y — > at x — > oo. 



5.2.1 Small baryon density 

We first study £(yo) in the regime where the baryon density is small. Since the full 
energy at finite /i takes the form 

£(yo,ti>) = £min(yo) - fiVo , 

it is important to determine the behavior of £ m i n (yo) for small y^. As long as the 
potential for /x = is quadratic (or linear) for small y^, we must have a first order 
transition to some finite baryon density at a critical chemical potential rather than 
a continuous transition where the baryon density increases gradually from zero. The 
results we obtain at small yo are also very robust (within our approximation), since here 
all field strengths and derivatives are small, and the incompletely known a' corrections 
in the D8-brane effective action are not important. 

The terms in (125]) coming from the Yang-Mills action are simply the leading order 
kinetic and potential terms, 

7 , 

CUq I f°° , ,1 5 / =7~w„ Jf„ \\2 , 1 \2 1 I- 




It is convenient to change variables to obtain a canonical kinetic term. Thus, we define 
u such that 

du 1 1 



dX xlyffix) ~ X " 

Choosing u = to correspond to x = 1, we have 

x(u) = sec3(-ii) . 

Note that x = oo corresponds to u — tt/3, so we now have a finite domain, which is 
convenient for our later numerical methods. Dropping the overall constant and working 
at /i = for now, we have 

£= 3 du{\{y'f + l -Xlx 2 {u)y i ) (26) 



o 



2 W/ 2 
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Extremizing, this gives rise to the differential equation 

y"(u) = 2\lx 2 {u)y\u) (27) 

As we discussed above, for a given y(0) > 0, solutions to this equation with slope 
larger or smaller than some critical value will approach positive or negative infinity 
as it — > 7r/3 and give rise to an infinite energy. The minimal energy configuration 
corresponds to the critical value of the initial slope for which the solution approaches 
zero at u = 7r/3. For y << 1/A, the solution is linear to a good approximation, since 
taking 

y(u) = y (l - -u) (28) 

7T 

we find that the right hand side of fl27|) is small enough that even the maximum value 
of y" integrated over the interval would only change y' slightly. 

Thus, for y << 1/A, the energy is given by inserting ( 1281) into ( 1261) . and we find 

3 

£ e ff{yo) ~ ttvI + 0{\ 2 Vq) small y 

Thus, the full energy S(y , fi) is always positive for small enough y , and the transition 
to nuclear matter must be first order in our model. 

While we can no longer trust the Yang-Mills approximation for large yo (of order 
l/y/X or larger), it still interesting to look at behavior of the Yang-Mills terms in the 
energy functional in this regime. Continuing to use only the terms (126]) . a numerical 
study suggests the asymptotic behavior 

1 „ 
4// ~ 3A02/0 

Note that this asymptotic growth in the energy density as a function of yo is not enough 
to stabilize the baryon density to finite values for jl larger than value 

A = ~A ■ 

Comparing with ( 1241) . we see that this value corresponds precisely to \i — M^. Thus, we 
conclude that the a! corrections in the Born-Infeld action are essential for stabilizing 
the baryon density to finite values for large /z, and that without these, the baryon 
density would diverge beyond a critical chemical potential that exactly coincides with 
the large A result for the baryon mass. In fact, we will see that at large A the Born-Infeld 
corrections only modify this critical chemical potential by terms of order -r. 

5.2.2 The critical chemical potential 

Now that we have demonstrated that there must be a first order phase transition to 
nuclear matter in our model, we would like to determine the critical value of /1 above 
which a non-zero baryon density is favored, and the baryon density as a function of fx 
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above this. Thus, we repeat our numerical study from the previous section, but this 
time with the full energy functional. In this case, the differential equation for y (using 
the same coordinates) is 

v " = -^ y ' f Tu + {{y ' )2 + ^ ln H ~ y Tu {l + HH) 

where 

H(x, y) = ^1 + A 2 (^ + 1(j,3_j,3 )2) . 

As before, the energy is minimized for a critical solution to this equation that ap- 
proaches at u = 7r/3. 

Our results indicate that the energy £ m i n (yo) behaves as a quartic function of yo for 
large yo, so the Born-Infeld terms stabilize the baryon density to finite values for any 
value of [A. As we have discussed, the critical value of the chemical potential beyond 
which a nuclear matter phase is favored is given by the minimum value of the energy 
per baryon. Specifically, we have 

£ 

/i c = mm yo ^; 

Vo 

We have numerically evaluated this critical chemical potential for large values of A 
ranging from A = 10 to A = 3000. Our data for fi crit at large A are fit very well with a 
function 

/W = M° B (A + cAq- 1 + 0{\~ 5 )) (29) 
where the best fit values are 

A « 1.000 ± 0.001 c^8±2 

Thus, to very good accuracy, the critical value of the chemical potential approaches 
the baryon mass for large A. Though our analysis using singular homogeneous configu- 
rations is an approximation, it is implausible that the almost exact agreement between 
the critical chemical potential and the baryon mass that we find here for large A is a 
numerical coincidence. A more plausible explanation is that the ratio of the critical 
chemical potential to the baryon mass does approach 1 in the limit of large lambda, 
and that our approximation gets this leading result correct. This is in accord with the 
expectation that our approximation should become exact in the limit of large baryon 
density, since as we will see below, the baryon density just above the transition does 
approach infinity as A becomes large. 

Thus, we believe that a robust conclusion of our analysis is that the binding energy 
per nucleon for large A is a vanishing fraction of the baryon mass. 

5.2.3 The binding energy per nucleon 

To determine the actual value of the binding energy, we need to compare the subleading 
term in f[2"9~l) with the subleading term in the baryon mass. 
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Even if our approximation is also correct for this subleading term, evaluating the 
coefficient c here depends crucially on the higher order terms in the Born-Infeld action. 
Since we have used the abelian D8-brane action together with an ad-hoc ordering 
prescription in lieu of the unknown full result for the effective action, we expect that 
the numerical value here is not reliable, However, the result that the correction is of 
order A -1 (rather than e.g. A~a) should be robust. 

Similarly, a correct calculation of d in the result 

M B = M B (l + ^ + ...) 
Ao 

for the baryon mass discussed in section 3.1 probably requires more complete knowledge 
of the non-abelian effective action. However, recalling that the leading order result for 
the baryon mass is proportional to A, we see that the result for the binding energy per 
nucleon (Mb — // c ) is actually relatively insensitive to A for large A. Since we also know 
that this binding energy approaches some constant value in the limit of small A (the 
large N c QCD result with two massless flavors), then assuming a smooth behavior at 
intermediate values of A, we can treat the large A result as a prediction for the order 
of mag nitude of the QCD resultQ 

Noting that Mkk ~ ^qcd for large A, the value of the binding energy per nucleon 
extrapolated to iV c = 3 becomes 

Eund = ^-^qcd{c' - c) « 7 MeV(c' - c) 
y7T 

As we have noted in the introduction, this is indeed of the same order of magnitude as 
the physical QCD result of lQMeV assuming that d — c is of order one. 

5.3 Baryon density above the transition 

We can also calculate the baryon density just above the transition. Our results suggest 
that just above the transition, the preferred value of yo for large A behaves like 

for K « 0.31. This suggests that 

n B R oc Aa 

as A is increased. This is consistent with the finding of Sakai and Sugimoto that the 
baryon size goes like A~a. 

For large chemical potential, the result that the \x = energy density approaches 

i 

t/o f° r large y implies that the baryon density minimizing S(y oc n B ) — jj,nB for large 
/i is 

ub oc fi 3 

16 Another example with similar insensitivity to A for both large and small A is the free energy of 
JV = 4 SUSY Yang-Mills theory. Here, it is indeed the case that the large A result for the free energy 
gives a good prediction of the order of magnitude of the the small lambda result (or vice versa) . 
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Also, the energy density as a function of \i for large \i behaves as 

E oc /i 4 

Note that the powers here are those appropriate for free fermions. We would like to 
understand this point better. 
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A Holographic Dictionary 

Consider a configuration (before decoupling) with a D4-brane in the 01234 directions 
(with £4 noncompact) and a D8-brane in the 012356789 directions, but bent in a U 
shape so as to intersect the D4-brane at two places along x 4 . Locally, one of these 
intersections is a D8 and the other is a D8. Now, we are interested in the coupling 
between the D8-brane gauge field and the operators 

B L = 

at the one intersection and 

Br = ip R ipR 

at the other intersection. We define i[) L and ip R such that they create particles with 
positive charge on the D4-brane, or physically, such that a test charge on the D4-brane 
is repelled from both of these particles. In this case, the baryon number operator is 

B = 1/N C (^ L + ^ R ) 

and we have a coupling a B in the effective action where ao is the time component of 
the D4-brane gauge field. 

We would now like to understand how the two operators and Br couple to 
the D8 brane gauge field. To do this, we note that if we perform a rotation by n in 
the 1-4 directions, centered at the point on the D4-brane between the two D8-branes, 
we get back to precisely the same configuration, since the D4-brane does not change 
orientation, while the D8 and D8 branes will switch orientation but also switch position. 

Now, suppose we have a configuration with one tjr L particle at the D8 intersection. 
This repels a test charge on the D4-brane, so after the rotation it is still a particle 
that repels a test charge, but now it is a particle at the D8-bar intersection. It must 
therefore be a ip R particle. Thus, a particle is mapped to ip R particle. Now, suppose 
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that we have a test charge on the D8-brane that is repelled by the particle in the initial 
configuration. In the rotated configuration, this test charge will still be repelled (by 
the ifr R particle). Also, the test charge in the new configuration will have the same sign 
as in the old configuration, since we have simply performed a rotation. This means 
that if we describe the entire U-shaped D8-brane using a single patch, both ijr L and %jr R 
particle will source electric fields pointing away from the D4-brane (or both towards 
the D4-brane, depending on our convention). 

This implies further that if we use a single field over the entire D8-brane config- 
uration, then the coupling of Aq to Bl at the D8 intersection will have the same sign 
as the coupling of A to Br at the D8-bar intersection. For the Sakai-Sugimoto setup, 
this implies that if we want to turn on a chemical potential for baryon number (i.e. 
turn on the operator B), we want to choose A to have the same sign at o = oo as at 
cr = — oo (if we use the same gauge field over the whole brane configuration). 
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